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Abstract 

Assuming three-neutrino mixing, we study the capabilities of very long baseline neutrino 



Oh 



oscillation experiments to verify and test the MSW effect and to measure the lepton mixing 



angle 613. We suppose that intense neutrino and antineutrino beams will become available 
in so-called neutrino factories. We find that the most promising and statistically significant 
results can be obtained by studying v e — > and v e — > oscillations which lead to matter 
^ ■ enhancements and suppressions of wrong sign muon rates. We show the #13 ranges where 
matter effects could be observed as a function of the baseline. We discuss the scaling of rates, 
significances and sensitivities with the relevant mixing angles and experimental parameters. 
Our analysis includes fluxes, event rates and statistical aspects so that the conclusions should 
be useful for the planning of experimental setups. We discuss the subleading Am^ effects 
in the case of the LMA MSW solution of the solar problem, showing that they are small 
for L > 7000 km. For shorter baselines, Am^ effects can be relevant and their dependence 
on L offers a further handle for the determination of the CP- violation phase 5. Finally we 
comment on the possibility to measure the specific distortion of the energy spectrum due 
to the MSW effect. 
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1 Introduction 



The long term aim to build muon colliders offers the very attractive intermediate possibility 
for "neutrino factories" Jl|, |2|, [!| with uniquely intense and precisely characterized neutrino 
and antineutrino beams. This requires only one muon beam at intermediate energies such 
that neutrino factories are rather realistic medium term projects which constitute also a 
useful step in accelerator technology towards a muon collider. The current knowledge of 
neutrino masses and mixing implies for typical setups very promising very long baseline 
neutrino experiments. We study in this paper in a three neutrino framework the potential 
to verify and test the MSW effect and to measure or limit #13 in terrestrial very long 
baseline experiments with neutrino factories where the neutrino spectrum and fluxes are 
rather well known and under control |T| . Calculating the oscillation probabilities and event 
rates for different channels and comparing with those for oscillations in vacuum we find 
that the asymmetry between the v e <-> and v e <-» oscillations is a very promising tool 
to test and verify the MSW effect. The reason is, as we will see, that matter effects lead 
to measurably enhanced event rates in one of the two channels and to equally suppressed 
transitions in the other. The asymmetry between the event rates associated with these 
two channels would therefore be very sensitive to the MSW effect since the matter induced 
changes have opposite effect on the two rates thus amplifying the "signal", while at the 
same time common backgrounds would drop out. We analyze event rates and we include 
statistical aspects such that the results are directly applicable for the planning of optimal 
experimental setups. We discuss the capabilities of a neutrino-factory experiment as a 
function of the distance between the neutrino source and the detector and of the muon 
source energy for the optimal observation of the MSW effect. Moreover, we determine the 
sensitivity to the value of #13 for different experimental configurations. 

Demonstrating and testing the MSW effect directly is of fundamental importance since this 
effect plays a basic role in different neutrino physics scenarios. The MSW mechanism pro- 
vides, for example, the only clue for understanding the solar neutrino deficit with a neutrino 
mass squared difference within a few orders of magnitude from that inferred from the at- 
mospheric neutrino data. Atmospheric neutrinos can undergo matter-enhanced transitions 
in the earth. The matter effects in neutrino oscillations will play an important role in the 
interpretation of the results of a neutrino factory experiment using a L > 1000 km baseline. 
They also play a role in the searches for CP- violation in such experiments 0, |], since 
matter effects generate an asymmetry between the two relevant CP-conjugated appearance 
channels (4| . Knowing the asymmetry caused by matter effects is therefore essential for ob- 
taining information on the CP-violation originating from the lepton mixing matrix. Matter 
enhanced neutrino transitions can play important role in astrophysics as well. 

Neutrino factories are extensively discussed in the literature |L], p|, |3|, m [7], §. Either 
muons or anti-muons are accelerated to an energy and decay then in straight sections of 
a storage ring like \T — > e~ + u e + or — > e + + + v e so that a very pure neutrino beam 
containing v e and or v e and z/ M , respectively, is produced. The muon energy could be 
in a wide range from 10 GeV to 50 GeV or more and a neutrino flux corresponding to 2 • 10 20 
muon decays per year in the straight section of the ring pointing to a remote detector could 
be achieved. Higher fluxes are also currently under discussion [0]. The neutrino fluxes are 
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therefore very intense and can be easily calculated from the decay spectrum at rest. For 
unpolarized muons and negligible beam divergence one finds for a baseline of L = 730 km 
a neutrino flux of ~ 4.3 • 10 12 yr _1 m -2 and for a baseline L = 10000 km a neutrino flux 
~ 2.3 • 10 10 yr _1 m~ 2 . Note also that the v e and v e fluxes depend sizably on the beam 
polarization, that will be assumed to be negligible in this paper. Altogether a neutrino 
factory would provide pure and high-intensity neutrino beams with a well known energy 
spectrum that in turn would allow a wide physical program including precise measurement 
of mixing parameters 0, ||, matter effects 0, § and, in case of LMA solution of the solar 
problem, leptonic CP- violation || [7j. 

The produced neutrino beam will be directed towards a remote detector at a given Nadir 
angle h, which corresponds to an oscillation baseline L = 2i? e cos(h) , where i2© = 6371 km 
is the earth radius. If L <; 10 3 km, as we shall assume, matter effects become important in 
neutrino oscillations^. For L < 10600 km the beam traverses the earth along a trajectory in 
the earth mantle without crossing the earth core where the density is substantially higher. 



According to the earth models |T^, [13[, the average matter density along the neutrino 
trajectories with L = (10 3 — 10 4 ) km lies in the interval ~ (2.9 - 4.8) g/cm 3 . The matter 
density changes along each trajectory, but the variation is relatively small - by about 1.5 to 
2.0 g/cm 3 , and even the largest takes place over relatively big distances of several thousand 
kilometers. As a consequence, one can approximate the earth mantle density profile by a 
constant average density distribution. The constant density should be chosen to be equal to 
the average density along every trajectory |TJL [TR [HI. For the calculation of the neutrino 



oscillation probabilities in the case of interest, the indicated constant density model provides 
a very good approximation to the somewhat more complicated density structure of the earth 
mantle (for a recent discussion see, e.g., [fT/jj). Let us note also that the earth mantle is 
with a good precision isotopically symmetric: Y e = 0.494 [T2|, 13], where Y e is the electron 
fraction number in the mantle. 

Since the beam consists always either of v e together with or v e in combination with z/ M 
there are, in principle, eight different appearance experiments and four different disappear- 
ance experiments which could be performed^. From an experimental point of view, however, 
at present the four channels with muon neutrinos or antineutrinos in the final state, — > 
Vfj, — > Pfj,, v e — > ^fi, v e — > Vfj, seem to be most promising. A very important issue is in this con- 
text the ability of the detector to discriminate between neutrinos and antineutrinos, namely 
the ability to measure the charge of the leptons produced by the neutrino charged current 
interactions. Note that very good discrimination capability is required for a measurement of 
the appearance probabilities since they produce "wrong sign" muon events in the detector 
and have to be discriminated from the much larger number of events associated with the 
and survival probabilities. The channels with electron neutrinos or antineutrinos in the 
final state are problematic from this point of view due to the difficulty of telling e + from e~ 
in a large high-density detector. On the contrary, a very good fJ. + /^~ discrimination could 
be obtained in a large properly oriented magnetized detector III 



1 The possibility to detect matter effects in long baseline neutrino oscillation experiments with L ~ 730 km 
(MINOS, CERN - GS) was discussed, e.g., in refs. @ 0. 

2 For pre-neutrino factory discussions of very long baseline terrestrial neutrino oscillation experiments 
see, e.g., |J). 
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The paper is organized as follows. In section 2 we give the analytic formulae for three 
neutrino oscillations in matter which contain the essential physics relevant for our study. 
In the following section 3 we discuss event rates, their parameter dependence, their scaling 
behaviour and we show results from our numerical calculations. In section 4 we define the 
sensitivity to matter effects in a statistical sense and discuss the results of our numerical 
calculations including parameter uncertainties. This is followed by a discussion in section 5 
of the effects of a non- vanishing Am^. In section 6 the possibility of detecting matter 
effects by looking for the characteristic enhancement, the broadening and the shift of the 
dominating maximum of the event rate energy spectrum is discussed and we conclude in 
section 7. 

2 Three Neutrino Oscillation Probabilities in Matter 

We will assume in this paper the existence of three flavour neutrino mixing: 

3 

Wi> =}Uik\vk > , l = e,n,r, (1) 
fc=i 

where \ui > is the state vector of the (left-handed) flavour neutrino ui, > is the state 
vector of a neutrino vj. possessing a definite mass m^, m^ ^ rrij, k ^ j = 1, 2, 3 and U is a 
3x3 unitary matrix - the lepton mixing matrix. We conventionally order the masses in such 
a way that < Am^ < [Am 2 -,]. According to the sign of Am^, the previous convention 
corresponds to the two closest neutrino masses (mi, 777.2) being lighter than the third one 
(777.3) (W3 > 777.2 > mi and Amjj > 0) or being heavier than the third one (7712 > mi > 777.3 
and Am^ < 0)Q. As long as matter and CP-violation effects are not taken into account, the 
two cases are phenomenologically equivalent from the point of view of neutrino oscillations. 
It is natural to suppose in this case that one of the two independent neutrino mass-squared 
differences, say Am^, is relevant for the vacuum oscillation (VO), small or large mixing 
angle MSW solutions^ (SMA MSW and LMA MSW) of the solar neutrino problem with 
values in the intervals [^(], |2lJ 

VO : 5.0 x KT 11 eV 2 £ Am 2 21 £ 5.0 x lO" 10 eV 2 , (2a) 
SMA MSW: 4.0 x 10" 6 eV 2 £ Am 2 21 £ 9.0 x HT 6 eV 2 , (2b) 
LMA MSW: 2.0 x 10" 5 eV 2 £ Am 2 21 £ 2.0 x HT 4 eV 2 , (2c) 

while Am^ is responsible for the dominant atmospheric <-> v T oscillations and lies in the 
interval 

ATM : lO^eV 2 £ (Am 2 ^ £ 8.0 x 10~ 2 eV 2 . (3) 

3 This scheme implies an approximate degeneracy of either mi and TO2 or of all the mass eigcnstates. 

4 There exists also a large mixing angle solution of the solar neutrino problem with Am^ ~ 10~ 7 eV 2 - 
the LOW solution, which, however, provides a somewhat worse description of the data than the MSW and 
the VO solutions ||]|. 
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For E v > 1 GeV, L < 10 4 km, the Am 2 -hierarchy 



Am 21 <C | Am 



31 1 



(4) 



and Amjj <C 10 4 eV 2 , the probabilities of three-neutrino oscillations in vacuum of interest 
reduce effectively to two-neutrino vacuum oscillation probabilities [^2, [23|]: 



COS 



Am^L 
2£ 



, l^l' = e,fi,T, (5) 



^ 3 : c (^-^)=^ c ^-^) = l-2|%| 2 (l-|%r) f 1 " 008 ^^). i = e,A*,r. (6) 

Under the conditions ([!]) and (|]) the element |Z7 e 3| of the lepton mixing matrix, which 
controls the v e — > i^( T ), ^ e ~~ ► ^mW' ~ ^ and the — > z/ e oscillations, is tightly 
constrained by the CHOOZ experiment and the oscillation interpretation of the solar 
and atmospheric neutrino data: for 3.0 x 10~ 3 eV 2 < Am^ < 8.0 x 10~ 3 eV 2 one has 

|l/" e3 | 2 £ 0.025 . (7) 

The CHOOZ upper limit is less stringent for 1.0 x 10" 3 eV 2 < \Am 2 31 \ < 3.0 x 10~ 3 eV 2 
where values of \U e 3\ 2 = 0.05 are allowed. 



Note that under the condition (Q), the VO or MSW transitions of solar neutrinos depend 

J 2 ,! and 6*i2, 



on \U e s\ 2 and on the two- neutrino VO or MSW transition probability with Am^ and Q\ 



where 



2o/. I l ( ' U < '-2\ 0/J _ \Uel\ 2 ~ \U £ 2 
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u e2 \ 


to 
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\u e2 


2) 2 



12 



sin 2 2g 12 = 4 ' 7' ■ ' P- - , cos2# 12 



|[/el| 2 + |t/e2| 2 



playing the role of the corresponding two-neutrino oscillation parameters |2"5] , |2*fi| , |2"Tf (for 
recent reviews see, e.g., []28], ^9|). For |?7 e 3| 2 satisfying the limit (^), however, this dependence 
is rather weak and cannot be used to further constrain or determine |t/ e 3| 2 . In general, under 
the condition (^) and for Am 21 <C 10~ 4 eV 2 the relevant solar neutrino transition probability 
depends only on the absolute values of the elements of the first row of the lepton mixing 
matrix, i.e., on |£7 e i| 2 , i=l,2,3, while the oscillations of the (atmospheric) v e and v e 

on earth distances are controlled by the elements of the third column of U, \Uis\ 2 , I — e,fi, t. 
The other elements of U are not accessible to direct experimental determination. Moreover, 
the CP— and T— violation effects in the oscillations of neutrinos are negligible. 

For our analysis we use a standard parametrization of the lepton mixing matrix U : 

— iS 

C12C13 S12C13 $136 

— SnC23 — Ci2S23Sl3e C12C23 — Si2S23 s 13 e S23C13 

(9) 

S12S23 - c 12 c 2 3Sne lS -C12S23 - si2C2 3 si 3 e J<5 c 23 ci 3 

where = cos%, s^- = sin 6^ and 5 is the Dirac CP- violation phaseQ. The angles 6\2 and 
#23 in (^]) are constrained to lie within rather narrow intervals by the solar and atmospheric 

5 We have not written explicitly the possible Majorana CP-violation phases which do not enter into the 
expressions for the oscillation probabilities of interest po| . Furthermore we assume throughout this study 
< 0i2, 023, 013 < tt/2, < S < 2tt. 
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neutrino data for each of the different solutions, eqs. (fja) - (0c), of the solar neutrino 
problem. With the accumulation of data the uncertainties in the knowledge of 612 and #23 
will diminish, while only upper limits on s 2 3 like eq. (0) have been obtained so far. The 
mixing angle #13 is one of the 4 (or 6 - depending on whether the massive neutrinos are of 
Dirac or Majorana type p0| , |23|j ) fundamental parameters in the lepton mixing matrix. It 
controls the probabilities of the v a {y^) — ► ^e(^), Vpipe) — > ^e(^v), — ► v T and v e — > u T 
oscillations and the CP— and T— violation effects in neutrino oscillations depends on it. 
Obviously, one of the main goals of the future neutrino oscillation experiments should be to 
determine the value of #13 or to obtain a more stringent experimental upper limit than the 
existing one @. 

Under the condition (f|), with Am^ <C 10~ 4 eV 2 and in constant density approximation, 



the oscillation probabilities of interest take the following simple form (see, e.g., p6l 51, [2 



Pf> e - v,) = s 2 23 Pf (Am^,sin 2 20 13 ) , 
K - ^) = 4 + 4 [1 - ^(Am^, sin 2 26 13 ) 



+ 2c 23 s 2 3 



cos k — (1 + cos 26*^) sin 



AE m L 



sin(K + 



AP m L, 



2sin 2 |-ip^(Am 2 1 ,sin 2 2£ 13 ) 



+ (1 + cos20 13 ) sin — - — sm(/«H — t 



where 



Pj^Am^, sin 2 26 13 ) = - [1 - cos AE m L] sin 2 20™ 



is the two-neutrino transition probability in matter with constant density and 

L 



K 



Am,, Tr . 

-if + y - AE * 



'10a) 



(10b) 



(10c) 



(11) 



(12) 



is a phase. In eqs. ( [0] ) - (|T2|), AP m and 6^ are the neutrino energy difference and mixing 
angle in matter, 



AP„ 



I Am 2 ! I 
2P 



1 / 2EV 
C + , =^ («■*■, -3^ 



(13) 



where 



and 



i=f 



2EV 



±4 



2EV 
Ami! 



sin 2 fl 



13 



1/ = v^Gf^ 



man 



(14) 



(15) 
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is the matter term, 



jvrman 

e 



being the average electron number density along the neutrino 



trajectory in the earth mantle, The probability P- 



3v, 



can be obtained from eq. (|T0| 



by replacing the factor s| 3 with c| 3 . The corresponding antineutrino transition and survival 
probabilities have the same form and can formally be obtained from eqs. ([10|a— |10]c) by 
changing the sign of the matter term, i.e. N™ an —>■ _j\r™ an in the expressions for k, AE m 
cos 26^, eqs. (|T2f)— (|13D , and by replacing C + by C_. Expressions (|T0|) follow directly from 
eqs. (12), (13b), (14b) and eq. (19) in |26[] (after setting in the corresponding formulae 
(p' 23 = 0, (p' 12 = vr/2). They can be deduced also from the explicit expressions for the 
probabilities of neutrino oscillations of the earth-core-crossing neutrinos in the two-layer 
approximation for the earth density distribution (see, e.g. 
(on the basis of |B6fl) and in 103 



14, 15, 16, Bl , obtained in m 



As we shall see in Section 5, at distances L < 6000 km and for Am^ ~ 10 4 eV 2 , the Am 2 , 



3i/ / 



corrections in the probabilities of interest P§ 



P 



3f I 



and Pl v (v. 



21 

can 



be non-negligible. Utilizing the results of it is not difficult to derive also the expressions 
for these probabilities including the leading order (CP-conserving and CP-violating) Am^- 



3/'/ 



corrections [33]. One finds in the case of the probability P|' 



p^) s 2 3 [1 + cos 5 cot 6 23 sin 26' 23 ] P%{ Am 31 , 6 13 ) 

AE L 

+ cos 6' 12 sin 2^23 sin 26™ sin — 



'13 



X 



cos 5 ( sinfft 



AEmL , 



cos 29^ sin 



AE m L 



n . c- . « . /« AE m L, 

zsmo sin- sin — I 

2 v 2 2 ' 



(16) 



where [26, 34 



Am 31 = Am 31 



'12 



Am 



21 5 



(17) 



K 



L 
2 



Am 



31 



2E 



V - AE ri 



LAml, 

21 cos 20 

2E 



12 



and 



cos 6^ 



12 



Am^ci 2 Si2C 13 



^ Am|]Ci2£l2Cl3 

2EV + s? 3 A^l - Am 2 cos 2fl 12 2 ^ + 4s Am li 



(19) 



sin 26o 



Am^sis sin 26*i2 



23 



Am 31 C| 3 — Am 21 cos2#i2 



Amji 
Am 2 x 



Si3 sin20 



12 



(20) 



The expressions for Am 31 ,0i3), Ai? m and cos 26^ are given by eqs. (JTT]) and ( |T3"D 



in which Am 31 is replaced by Am 31 . The probability P^{y e 



(Pi 



3i/ / 



z/ e )) can 
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be obtain from eq. ([16]) by making the change V — > — V and 5 — > — 5 (5 — > — 5). 
The expression flTBD was obtained by neglecting the ~ (cos^) 2 , (sin20 23 ) 2 , cos 9' 12 sin 29' 23 
and the higher order corrections (for further details see [[26], [34]]). Note that if Am^ < 
2 x 10~ 4 eV 2 and for iV e man > 1.45 cm- 3 N A (L > 1000 km), we have 2EV > 10~ 3 eV 2 } and 
consequently | cos#' 12 | < 0.1 at E > 4.5 GeV. For the indicated maximal value of Am^ and 
lAmfil > 10~ 3 eV 2 , s 2 3 < 0.025 (0.05) one finds | sin2fl 23 | < 0.032 (0.045). It follows from 
the preceding discussion that if (Am^l = 3.5 x 10~ 3 eV 2 , the correction ~ cos9' 12 w °uld be 
the dominant one for E ~ (5 — 30) GeV. For the values of Am^ and Am^ of interest, the 
Am^- correction included in Am 31 , eq. (JT7|), essentially just shifts the resonance energy by 
at most ~ 10%. The terms with the factors cos S and sin<5 in eq. (16) include the leading 
order CP-conserving and CP- violating Am^- corrections, associated with the phase 5. 

Several comments are in order. First, the analytic expressions eqs. (|10]a)— (|10|c) represent 
excellent approximations in the case of the VO and SMA MSW solutions of the solar neutrino 
problem and for values of Am^ ^ 5 x 10~ 5 eV 2 from the LMA MSW solution region. For 
5 x 10~ 5 eV 2 < Am 21 ^2x 10~ 4 eV 2 , the corrections due to Am 21 can be non-negligible and 
we are going to discuss them in section |5|. Second, as it is clear from eqs. (|TU|a) and (|TTJ) , the 
probabilities P^[y e — > v^) and P^[p e — > z/ M ) cannot exceed s 23 ~ 0.5. The maximal value 
s 23 can be reached only if the MSW resonance condition sin 2 26^ = 1 and the condition 
cosAE m L = —1 are simultaneously fulfilled. At the MSW resonance, however, one has 
AE res ^ 12 3 7r x io~ 4 km" 1 tan 2#i 3 A e man [iV A cm- 3 ], where iV e man is in units of iV^cm -3 , N A 
being the Avogadro number, and for s 2 3 < 0.025 the second condition cosAE^ s L = —1 can 
only be satisfied for L > 10 4 km. If s 2 3 = 0.05 this condition requires L > 8 x 10 3 km. Given 
the fact that the neutrino fluxes decrease with the distance as L~ 2 , the above discussion 
suggests that the maximum of the event distribution in E and L in the transition v e — + 
should take place approximately at the MSW resonance energy, but at values of L smaller 
than those at which max{P^[y e — > u^)) = s 23 . The MSW resonance energy is given 
by: E ms = 6.56 x Am 31 [10~ 3 eV 2 ] cos 26'i 3 (A"™ an [cm -3 A^^]) -1 GeV, where Am 2 3l is in units 
of 10" 3 eV 2 . For Am^JlO" 3 eV 2 ] = 3.5; 6.0, 8.0 and, e.g., N™ an [cm- 3 N A } = 2, we have 
E res = 11.5; 19.7; 26.2 cos 2#i 3 GeV. The v e —> and v e — > oscillations will be affected 
substantially by the MSW effect if the energy of the parent muon beam E^ > E rcs . This 
condition can be satisfied for any value of Am 31 from the interval @ for E^ ^ 30 GeV. For 
the value of Am 31 [10 -3 eV 2 ] = 3.5 which is currently "preferred" by the Super-Kamiokande 
data we have ^ 15 GeV. Note that a change in the sign of Am 31 reverses the sign of 
N^ es which opens the possibility to determine this sign via matter effects, as was noticed 
recently also in ref. ||. The probabilities (|T0|) for L = 7330 km (the Fermilab - Gran Sasso 
distance) have been discussed recently in [|T7|]. 



3 Event Rates 



We discuss now the effects of matter on the total rate of u e — > v e — > P M , — > — > 
events. If, for example, fi + are accumulated in the storage ring, the neutrino beam contains 
v e and v^. Since the neutrino-antineutrino transitions have a negligible rate, the beam- 
induced /x + events in the detector must be attributed in this case to the charged current 
interactions of unoscillated v^. The total number n M + (/i + ) of ji + events measures therefore an 
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averaged survival probability. Wrong sign /j~ events must be attributed to generated 
by oscillations of the initial u e , so that their total number n M +(/i~) measures the averaged 
v e ~^ V/j, oscillation probability. Analogously, for fi~ in the storage ring, the total numbers 
of fj,~ and n + events ^-(/i - ), n^-(n + ) measure the averaged survival probability and 
—>■ oscillation probability, respectively. 

The total number of events in each channel are given by 

i q9/y e 3 r E » 

n M+ (/0 = N^N^—^-t / f v ^{E)Pl% e - u,){E) (dE/E,) , (21a) 
rv(/i+) = AT^iV kT l^S [*" f p ^(E)P 3 ^^^(E)(dE/E,) , (21b) 



n, 



in 9 Na e 3 r E » 

(A!") = iViV kT ^i-f / f^(E)P 3 ^ -> ^(E) (dE/E M ) , (21c) 
Tn ,,71 lj I v? . 



^ + (/i + ) = AT M _iV kT ^^^ / /^(E)Pf(* M -+ p„)(E) (dE/E,) , (21d) 



ft 



where Pg u denotes the oscillation probability described in section ||, iV M + (iV M - ) is the number 
of "useful" /i + decays, namely the number of decays occurring in the straight section 
of the storage ring pointing to the detector, iV kT is the size of the detector in kilotons, 
10 9 Na is the number of nucleons in a kiloton, E M is the energy of the muons in the ring and 
E m i n = 3 GeV is a lower cut on the neutrino energies that helps a good detection efficiency 
Since low energy events are suppressed by the low initial flux and the low cross section (see 
below), the results do not depend significantly on the precise value of E min for E min < 5 GeV. 
The functions / averaging the probabilities are given by 

f v&Vfl (E) = g Ve (E/E M ) [a Vll (E) /E^)e^- (E) , (22a) 

U^{E) = g D XE/E,){a^{E)/E^ + {E) , (22b) 

= g^E/E,){a^{E)/E,)e,-{E) , (22c) 

U^{E) = g^(E/E,)(a^(E)/E,)e, + (E) (22d) 

and take into account the appropriately normalized initial spectrum of i/-neutrinos produced 
in the decay of unpolarized muons, g v (E/E^), the charged current cross section per nucleon, 
<t v jp ){E), and the efficiency for the detection of \i~ (/i + ), e^-u+^E) (we neglect here the 
finite resolution of the detector). For the numerical calculations we use 

9vX x ) = 9vX x ) = 12x 2 (l - x) , g^(x) = g Dfi (x) = 2x 2 (3 - 2x) , (23a) 

a Vyk (E) = 0.67 • 10" 38 E cm 2 /GeV , (E) = 0.34 • 10~ 38 E cm 2 /GeV (23b) 

and e M -(E) = e M + (E) = e for E > E min so that /„„„„(£)/ U e ^(E) = f^(E)/f p ^(E) = 2 
independently of the energy. 

The contribution of the background to the number of muons observed in the detector has 
been neglected in eqs. ([H]). That background includes muons from the decay of charm 
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quarks produced by charged and neutral current neutrino interactions in the detector and 
from the decay of r produced by v T interactions. Both these sources can be kept under 
control in different ways |35], Qj. Let us note also that in eqs. (|2l| ) the possible divergence 
of the muon beam in the straight section of the storage ring has been neglected. 

We can analyze now the dependence of n^+([i~), n^-(fi + ), n^- (//"), on the mixing 

angles, the CP-phase 8 and some of the experimental parameters. In the Am 2 2 *C 10~ 4 eV 2 
approximation discussed above, the only relevant mixing angles are 823 and #13. In the 
leading order in a #13 expansion the dependence of the total event rates in matter (and in 
vacuum) on these parameters is 

/^-(/OjfV+C/^) ^ sin 2 26>23 , (24a) 

ra M +( J tr),n jU -Gu + ) cx sin 2 ^3 sin 2 2^13 . (24b) 

Higher order corrections in # 13 are constrained by the CHOOZ limit. Despite the resonant 
matter enhancement of the mixing due to #13, such corrections become only sizable for 13 
close to its upper limit and very long baselines where they can reach about 20% in the 
resonance channels. 



Since eq. (p4]b) will turn out to be useful when discussing the sensitivity to matter effects, 
we discuss it in greater detail. Let us write eqs. (|l0|a) in the form 

Pf> e - u„ v e -> P M ) = sin 2 8 23 sin 2 26 13 " , (25) 

where C+ is given by eq. ([[4]) and corresponds to neutrinos, C_ to antineutrinos and A 3i = 
&m\ x L/{AE). Note that C + enhances the rates when Am 2 x > 0, depletes them when 
Am^ < 0. Therefore the sign of Am 2 x determines whether the enhanced channel is the 
neutrino or the antineutrino one. In the limit in which sin 2 # 13 can be neglected on the 
right-hand side of eq. (|I4]) , eq. ( p5|) shows that ?v(aO> n^-{ji + ) are indeed proportional 



relevant when the first term vanishes around the resonance. For 2EV = Am| x we have in 



to sin 26* 13 . Despite the CHOOZ limit, sin 9 13 < 0.025, the second term in eq. ( |14|) can be 
relevant when the first 
fact C + = 2 sin #13 and 

q„/ x •>„ 9 „ L 2 V 2 ( sin(sin 9 13 LV) . 

Pf> e - Vfk ) = sin 2 £23 sin 2 26 13 — - \ 13 } , 26 

4 \ sm&i 3 LV 



whereas by neglecting the sin #13 term in fll4[) we would get C + = and 

Ff{y e -> = sin 2 6 23 sin 2 26 13 ^ . (27) 

A comparison of eqs. (B^) and ( |2"T| ) shows however that this approximation works even at the 
resonance provided L < 7r/(4\^sin # 13 ) ~ 7000 km (0.15/ sin# 13 ). The effect of the sin 2 # 13 
term in eq. fll~4] ) is therefore maximal in eq. ( |i~4"|) for very long baselines and #13 close to 
the experimental limit. In this case it can affect the rates sizably, while it is negligible for 
smaller #13 or smaller baseline. For a better approximation one can use, for muon energies 
such that the main contribution to the rates comes from neutrinos at the resonance, 



> / sinfsin 9i 3 LV) \ 
/T ,7V /i + ocsin 2 fl23sin 2 2g 13 . a " - » 2 ^ 
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which deviates less than 5% from the exact result in the whole parameter space (L < 
10000 km, 20GeV < < 50GeV). In most cases L < ir/(4V sin 13 ) holds and the 
oscillating term in eq. (p8| ) can be expanded, giving eqs. (H) or higher order approximations. 
We stress that in eq. (|28|) the baseline L appears as part of the correction to the 6*i 3 , 6*23 
scaling only. The dependence on L of the rates is more involved, especially for very large 
baselines, and will be described in figs. [I], ||. Note that the dependence on the beam intensity, 
detector size and efficiency is trivial: 

, n^- oc iV M + N kT e M + , (29a) 
?V(/i + ),rV + (/i + ) oc iViV kT e M - • (29b) 

We present now quantitative results for the total rates in matter for Am^ <C 10~ 4 eV 2 and 
we compare them with the results one would obtain in vacuum. The statistical significance 
of the matter effects will be discussed in the following section and effects of larger Am^ 
will be covered in section ^. 

The total event rates depend, as already discussed, in a transparent way on the experimental 
parameters N^±, iV kT , e M ± and on the mixing parameters #23, #13. We focus our discussion 
therefore on the less transparent dependence on the baseline and muon energy. For that 
we use the central value of Am^ = 3.5 • 10 _3 eV 2 , and we assume Am^ > 0. In the 
Anijj <C 10 -4 eV 2 approximation, in which CP-violation effects are negligible, the results 
for Am^ < can be obtained by simply interchanging neutrinos and antineutrinos. 

The total number of events in the two appearance channels v e — > v^, v e — > is shown 
for Efj, = 20 GeV and E^ = 50 GeV in fig. [I] as a function of the baseline (solid lines) 
in comparison with the event rates one would get if the neutrinos did not interact with 
matter (dashed lines). Fig. ^| shows the same results for the two disappearance channels 
u n — ¥ *V — > Vfj,. Both figures correspond to a "default" experimental set-up providing 
Nfj,+ = Np- — Nfj, — 2 ■ 10 20 useful muon decays (e.g. in one year of running) and to a 
detector with iV^r = 10 kilotons and an efficiency e M + = e^- = e = 50% in both channels. 
The rates depend only on the combination N^N^e which is in our case N^N^e = 10 21 . 
We assume for these figures that sin 2# 23 = 1 and sin 2 2# 13 = 0.01, one order of magnitude 
below the experimental limit. The rates for different values of N^, iV k T, e, ^23; ^13 can be 
obtained by using eqs. (p9|) , (0a) and (p4|b) or (|28|) . 

The vacuum event rates in the neutrino channels are twice as big as the rates in the an- 
tineutrino channels. This is because the oscillation probabilities are in the CP-conjugated 
channels in the Am^ 10 _4 eV 2 approximation the same while the functions averaging 
the probabilities are larger by a factor 2 in the neutrino channels (due to the larger cross- 
section). The disappearance channels shown in fig. |2| are essentially independent of matter 
effects since these effects come only with the #13 corrections to the oscillation probabili- 
ties^. In contrast, fig. [I] shows the drastic enhancement (depletion) of the event rates in the 
[p e — * Ufj) channel for very long baselines. The growth of the total rates with the 
muon energy which is obvious from eqs. (|2l|) can also be seen in fig. [I]. 

6 As discussed above effects as iarge as 20% can occur for 813 close to the CHOOZ limit and very large 
baselines. 
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0.7 1 2 3 4 5 6 7 10 0.7 1 2 3 4 5 6 7 10 

L/(1000 km) L/(1000 km) 



Figure 1: Appearance event rates n M +(/i~), n^-(jj + ) due to the u e — > and u e — > 
transitions in the earth mantle ( solid lines ) and in vacuum ( dashed lines ), as functions of 
the baseline L for = 20 GeV (left) and = 50 GeV (right). The differences between 
the solid and dashed lines is a measure of the magnitude of the earth matter effect. Both 
plots assume — 2 ■ 10 20 , e = 50%, sin 2 2023 = 1 and sin 2 26*13 — 0.01. The scaling of the 
rates with these parameters is described in the text. 



E^ = 20 GeV E„ = 50 GeV 




0.7 1 2 3 4 5 6 7 10 0.7 1 2 3 4 5 6 7 10 

L/(1000 km) L/(1000 km) 



Figure 2: Same as in fig. [J but for the disappearance channels. The dashed lines coin- 
cide almost perfectly with the solid lines, showing that matter effects are negligible in these 
channels. 
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Figs. |3] and |] show in more detail the differential event rates of the appearance channels 
v e ~^ an d v e Up, as three dimensional plots over L and E v . The figures show nicely the 
enhancement (suppression) due to the MSW mechanism in the v e — > {p e — > v^) channel 
for large baseline L. To understand these figures we can look at the values of dn^+(fi~)/dE 
and dn ll ~(fi + )/dE for three different energies, E = (5; 10; 15) GeV, and two baselines L. 
Specifically we compare a relatively short baseline L = Lq = 1000 km with a long baseline 
L = 6000 km where matter effects are more important. For L = 1000 km (L = 6000 km) the 
average matter electron density is N™ an = 1.45 ciii" 3 Na (_/\r™ an ^ 2.0 cul~ 3 Na), resulting in 
a resonance neutrino energy of E res = 15.8 GeV (E res = 11.5 GeV). The resonance energy 
is thus only somewhat reduced for larger baselines. The corresponding values of A 3 i(L) = 



dn^+(p-)/dE (GeV -1 ) dn^-(p+)/dE (GeV -1 ) 




Figure 3: Differential appearance rates for the channels v e — > (left) and v e — > (right) 
as functions of L and E for the same values of the parameters for which fig. [J was ob- 
tained. The asymmetry, i.e. the enhancement of dn^+(fi~)/dE and the suppression of 
dn^-i/i^/dE, atL> (5000 - 6000)km is related to the MSW effect (see the text). 



dn fl+ (jj-)/dE (GeV -1 ) dn^-(p+)/dE (GeV -1 ) 




Figure 4: Same as in fig. [| but for a beam energy of 50 GeV. 
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L Am| 1 /(4£?) and the matter effect factors C + and C_ which enter in P^{y e —* z/ M ) and 
PsiPe -> ^) (see eqs. (0), (0), (0) and (H)) are displayed in Table |. 



L 


1000km 


6000km 


E 


5GeV 


lOGeV 


15GcV 


5GeV 


lOGeV 


15GeV 


Agi(L) 


0.282tt 


0.1417T 


0.094tt 


1.6947T 


0.847tt 


0.565tt 


^ + 


0.686 


0.371 


0.111 


0.568 


0.159 


0.327 




1.315 


1.630 


1.945 


1.434 


1.869 


2.304 


A 31 (L)C+ 


0.610 


0.167 


0.031 


3.023 


0.424 


0.580 


A 31 (L)C_ 


1.166 


0.723 


0.575 


7.632 


4.973 


4.090 



Table 1: Values of A 31 (L) = LAm 2 31 /(4:E), C + and C_ for Am 2 31 = 3.5 x 10- 3 eV 2 , 
sin 2 26 23 = 1 and sin 2 26 13 = 0.01 (see the text). 



In general, A 31 (L) decreases linearly with the increase of E. At 1000km matter effects are 
small. C + (eq. [14]) decreases as E changes from 5 to 15 GeV, approaching the resonance 
energy from below. Note that due to the matter term we have C + < 1 and C_ > 1 for the 
three values of E of interest. The argument of the oscillating sine factor in P^{y e — > v^) 
and in P^{p e — » v^) at values of E considered is relatively small allowing an expansion of 
the sine in the oscillation probability P^{y e — > z/ M ) = s 2 3 sin 2 (2#i 3 ) (A 31 (L )) 2 . The fact 
that the matter effect factors cancel out reflects the negligibility of matter effects at short 
baselines[|. Neglecting matter effects at short baselines, the maximum of the differential 
event rate spectrum is at one half of the muon beam energy which is 10 GeV in the case 
under discussion. 

At L — 6000 km, one has $(L) = $(L )(L /L) 2 = $(L )/36, where $(L) is the flux of u e 
or i> e at distance L > L$. The argument of the sine is large for the anti-neutrino channel 
giving for the corresponding rates the upper bound 

dn fl ~(fi + )/dE oc s 2 3 g Pe (x) sin 2 (A 31 (L)C_) 

2 $(Lq) sin 2 29 13 
< s 23 g Pe {x) — — 2 — . (30) 

Obviously, dn^- (/j, + )/dE is strongly suppressed primarily by the decreasing of the flux. The 
suppression due to the matter effect term (1/C_ 2 ) for E = (5 — 15) GeV is by a factor of 
~ (2 — 5). The same conclusion is valid for dn ll +(fi~)/dE below the resonance region, say 
at E = 5 GeV. 

For neutrinos in the resonance region E = (10 — 15) GeV, the argument of the sine is small 
even for very long baselines. Like in the case of short baselines an expansion becomes 



7 Note however that for longer baselines where A <C 1 does no longer hold, the same conclusion would 
be wrong. 
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possible: 



dn„ + (»-)/dE oc 4 g Ve (x) $(A>)(^) 2 sin 2 (A 31 (L)C + ) 



4 sin 2 (2^ 13 ) (-^ L ) 2 . (31) 

This is approximately equal to the same differential event rate at short baselines. Thus, in 
the resonance region {E = (10 — 15) GeV) and at sufficiently long baselines (L ^ 4000), 
matter effects lead to a strong suppression of anti-neutrino event rate while keeping the 
neutrino event rate essentially constant with the change of L from ~ 1000km to ~ 7000km 
(see fig. |H). In vacuum, however, the neutrino rate as well as the anti-neutrino rate would 
show the same strong (L /L) 2 suppression at long baselines. 

For fixed L > 4500 km, say at L = 6000 km, the matter effects lead to a strong enhancement 
of dn fJi +([i~)/dE as a function of E in a region which is somewhat wider than the MSW 
resonance region, E ~ (8 — 15) GeV, in spite of the fact that in this region A 3 i(L)C + < 0.6 
and 



dn^-)/dE oc 4 g Ue (x) ^ 2 sin 2 (A 3 i(L)C 



sin 2 26»i3 

c~ 2 

s 2 23 g Ue (x) $(L) sin 2 2^ (A 31 (L)) 2 . (32) 



We have in the indicated energy region A31 (L) > 1 and the enhancement of dn^+ (fx ) / dE 
compared to the case of v e «-> vacuum-oscillations is actually given by the ratio of 
(A 31 (L)) 2 andsin 2 (A 31 (L)). 

Fig. [| shows the results corresponding to fig. |3] for a beam energy of 50 GeV. At very long 
baselines, the neutrino rate spectrum clearly peaks at the energy where the MSW-resonance 
condition is fulfilled. The change of the shape of the spectrum with the baseline L is more 
pronounced now because the resonance energy does not coincide with the maximum of the 
beam spectrum as it approximately does for a beam energy of 20 GeV. The baseline and 
muon energy dependence of matter effects will be further discussed in the next section in 
connection with a quantitative analysis of the significance of the effects shown in figs. [I] and 



4 Statistical Significance of Matter Effects 

We have seen in the previous section that matter effects change the total event rates in the 
appearance channels v e — > and v e —>■ in very long baseline experiments in a drastic 
way. Such experiments would therefore offer unique possibilities to observe matter effects 
and to test the predictions of the MSW theory. In order to study the capabilities of a 
neutrino factory experiment quantitatively, we must first define the meaning of "observing 
matter effects" . One of the most interesting possibilities would be a detailed observation of 
the shape of the neutrino energy spectrum which is modified by the MSW effect in a very 
characteristic way. This would allow to test non-trivial predictions of the MSW theory and 
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would allow to unambiguously attribute the enhancement/depletion of the total number of 
neutrino events to matter effects. High differential event rates and a good calibration of 
the detector would however be necessary for this option. We will discuss this possibility 
in sect. ||. In this section we confine ourselves to a discussion of the significance of MSW 
effects in total event rates. 

Matter effects produce deviations of the total number of wrong sign muon events n fl +(fi~), 
n>fi- from what is expected in the absence of matter. The discussion above clearly shows 
that such deviations occur in opposite directions in the appearance channels. Suppose that 
a certain number of wrong sign muon events in the neutrino (antineutrino) appearance 
channel n A1 +( / u _ ) (n At -( / u + )) is measured, while (n^+ c (/i~)) ((n™ c would be expected in 
the absence of matter effects. We want to determine the confidence level at which n fl +(fi~) 
and n^-(fi + ) could represent statistical fluctuations around the expected values in vacuum 



(n™ c (/i )), (n™ c (/j + )). We follow the procedure proposed by the Particle Data Book p6 | 



njAv ) 

(33) 

v- (^ + ) 



and calculate confidence levels by using 

X 2 = 2 [(n; a + c (/i-)> - n M+ 0O] + 2rv + (/0 log 

+2 [«_ c (/i + )> - v- (//+)] + 2rv log 



The corresponding "number of standard deviations" is given by n a = y x 2 - This prescrip- 
tion incorporates the available information in both the measured numbers n M +(/i~), n^- 
in the most complete way. Note that the results which will be shown later with this method 
in figures |5|-|7] below assume that #13 is known with some precision before such an analysis 
is performed. This could e.g. be realized by a measurement of the v e — ► and/or v e — ► z/ M 
rates at a relatively short baseline, say, below 1000 km. If #13 were only limited from above 
then one would have to vary #13 in the vacuum rates entering the definition of \ 2 m ec L- ( p3| ) 
in the allowed range and find the minimum. This would reduce the significance at largest 
baselines. The statistical method defined above in eq. (B3) can be understood intuitively 
by defining an asymmetry 

£ rv + (M-)+2rv(/x+) 1 ] 

which is obviously very sensitive to effects which affect n^+(/i~) and n At -(/i + ) in a opposite 
way, while at the same time a number of common systematic effects drop out. For small 
E vac (i.e. small beam energies) this method is related to comparing the absolute asymmetry 
expected in matter with the expected fluctuations of this quantity in the vacuum-case 
X = A/5A vac . For large E vac this is equivalent to doing the same with the asymmetry: 
X = A/5A vac . Note that the sign of Am^ can be determined from A or equivalently from 
A. This can easily be seen by observing, for example, that for = AL- the difference 
P E V ( u e -> v^-Pgipe -> u^, with Pf^e -> v,i) and P'f{pe -> Vy) given by eq. @, enters 
in A. The sign of A and A depends therefore unambiguously on whether C + or C_ is bigger 
than one (and the other smaller than one). Via the definition of C± given in eq. (|HD this 
is in turn unambiguously related to the sign of Am^. 
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Before we show the numerical results, let us discuss the qualitative dependence of n a on 
the relevant parameters. The dependence of n a on the intensity of the muon source and 
the detector size and efficiency, as well as the dependence on the mixing parameters follows 
simply from the previous section: 

n a (x (N^N^e^) 1 ! 2 and n a oc sin #23 sin 29\z . (35) 

The dependence of n a on the baseline L and the muon energy is less trivial. In gen- 
eral, matter effects increase with the baseline L. Thus a long baseline is essential for the 
observation of matter effects. For relatively small baselines L < 730 km, in fact, matter 
effects mostly cancel in eq. (^), thus making such baselines better suited for CP- violation 
measurements [0, [TJ . The dependence on needs more explanation and there is also some 
dependence on the statistical method used. In our simple approach which takes into account 
only total rates, the significance to matter effects grows first strongly with the increase of 
the beam energy. At very long baselines, where matter effects are non negligible, a maxi- 
mum is reached at a beam energy which allows a maximal number of neutrinos to fall into 
the MSW-resonance energy region. After a small depletion the significance starts to grow 
again. This is due to the fact that the rates of high energy (> 50 GeV) neutrinos and anti- 
neutrinos are strongly suppressed due to the matter effects. For large beam energies it is, 
as was already mentioned earlier, crucial to know the precise value of sin 2 26*13. Otherwise 
it is difficult to distinguish matter suppression from a reduced event rate due to a smaller 
6*i3 in absence of matter effects (see fig. [1] for large L). 

The dependence of n a on L and E^ is shown in fig. ^ for Am| x = 3.5 ■ 10~ 3 eV 2 and 
N^N^e = 10 21 , where the contour lines of n a are plotted in the L-E^ plane. The solid 




Figure 5: Contour lines of n a in the L-E^ plane. We use Am^ = 3.5 • 10~ 3 eV 2 and 
N^Nkje = 10 21 and the solid contour lines correspond to n a = 10sin26 ) i3 • {1,2,4,8,16}. 
The left plot assumes that #13 is known, while the right plot is obtained by varying #13 in the 
range currently allowed. 



16 



contour lines correspond to n a = 10 sin 2013 • {1,2,4,8, 16}. The left plot assumes that 6*13 
is known, while the right plot is obtained by varying # 13 in the range currently allowed. 
With the simple "total rates based" method described above and assuming in the following 
discussion that #13 is known, any L from the interval (4.0 — 10.0) x 10 3 km would be 
suitable for the purpose if Am 2 31 = 3.5 x 10" 3 eV 2 , sin 2 26 13 = 0.01 and = 20 GeV, 
with the sensitivity to matter effects increasing with L. For = 50 GeV, the sensitivity 
varies very little for L = (4.0 - 10.0) x 10 3 km. Let us note that at L > 3.0 x 10 3 
km the matter effects are substantial for any neutrino energy E either in the v e — > 
and/or in the v e — > channel, depending on the sign of Am^: at L > 3.0 x 10 3 km we 
have y/2G F N™ an L/2 > 0.95. The matter effects can be negligible for certain intervals of 
neutrino energies for L < 2.0 x 10 3 km. For relatively large E and L < 2.0 x 10 3 km, the 
oscillating factors in P^{y e — > v^) and P^{p e — > z/^), eqs. (pTOj) , (|TT|), (p~3| ) and ([14]), can 



be expanded in power series of their arguments and matter effects cancel in the leading 
orders of these expansions. Moreover, one has sin 2 (A 31 (L)) = (A 31 (L)) 2 . At L = 10 3 km 
and for Am^ = 3.5 x 10~ 3 eV 2 , sin 2 26*13 = 0.01, for instance, the earth matter effects 
cause a difference between P^{y e — > z/ M ), eq. flTU|), and the probability of the v e — > z/ M 
transitions in vacuum P 3 ^ c (z/ e — > i/ M ), eq. (^|), which does not exceed ~ 5 x 10~ 5 at E > 10 
GeV. We have, however, \P%{y e -> P M ) - Pf ac {p e -> z^)| < 10" 4 (10~ 5 ) for E > 10 (20) 
GeV. For Am 2 31 = 7.0 x 10~ 3 eV 2 , the same conclusions hold for approximately two times 
higher neutrino energies. Note that P^{p e ~^ u n) ~ S23 sin 2 26*13 an d for the values of 
sin 2 2#i3 and S23 used in the above examples we have s?, 3 sin 2 2#i3 = 0.005. Thus, baselines 
shorter than about 2000 km and parent muon energies E^ > (20 — 25) GeV seem to be 
better suited for searches of CP-violation in neutrino oscillations whose source is the lepton 
mixing matrix. The magnitude of the asymmetry between the probabilities due to matter 
effects effectively grows quadratically with the baseline until the argument A 31 (L)C + and/or 
A31 (L)C- approaches and exceeds one 0. 

The significance of the effect, i.e. the number of standard deviations n a) depends most 
crucially on #13 and L. We illustrate this dependence in fig. ^again for Am^ = 3.5-10~ 3 eV 2 , 
where the contour lines corresponding to n a = 1,2,3,4,5 are plotted in the sin 2 2#i3-L 
plane for two values of the product N^N kT e: 10 21 (left plot) and 10 22 (right plot). The 
muon energy is in both cases of fig. || E^ = 20 GeV, while fig. [7| shows the same plots 
with identical parameters for E^ = 50 GeV. The vertical dashed lines represent in all these 
figures the upper limit on sin 2 26*13. Figs. § and [7|show that matter effects could be observed 
in the total event rates for given baseline L in a rather large sin 2 2#i3 interval, while non- 
observation implies very strong upper bounds on sin 2 2# 13 . Figs. |6| and [7] show in other 
words the sin 2 26 l 13 range where the enhancement/depletion of the total appearance rates in 
vacuum due to matter effects is statistically significant for a given confidence level. In those 
ranges one can not only observe the deviations from the results expected in vacuum, but 
also the deviation from the results which one would get in matter if the sign of Am| x were 
reverted. A measurement of the sign of Am| x would therefore be possible in the sin 2 2# 13 
range where matter effects are statistically significant at a given confidence level. 

Finally we show in figs. |8| and ^ the sensitivity of the statistical significance to the chosen 
value of Affljj. Fig. |^ is exactly the same plot as fig. |5] with minimal and maximal Am 31 . 
In the left plot we have Am 2 x = 1.0 • 10~ 3 eV 2 and in the right plot Am 2 x = 8.0 • 10~ 3 eV 2 . 
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Figure 6: Contour lines of n a corresponding to n a = 1,2,3,4,5 in the sin 2 20i 3 -L plane 
for = 20 GeV and the two different values N^N^e = 10 21 (left plot) and N^N^e = 10 22 
(right plot). 




Figure 7: Same as in fig. [6| but for = 50 GeV. 
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Am 2 = 1.0 1(T 3 eV 2 Am 2 = 8.0 10" 3 eV 2 
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Figure 8: Sensitivity of the statistical significance of matter effects to the value of Am| x 
analogous to fig. [|. The left plot uses Am^ = 1.0 • 10~ 3 eV 2 while for the right plot Am^ = 
8.0 • 10~ 3 eV 2 with otherwise unchanged parameters. As in fig. [| the solid contour lines 
correspond to n a = 10 sin 2#i3 • {1,2, 4, 8, 16}. 




Figure 9: Sensitivity of the statistical significance of matter effects to the value of AmJ 
for fixed L = 8000 km as a function of E^. The lines show 10 sin 2#i3n CT for different Am| 
values. 
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Fig. |9] shows the statistical significance (i.e. 10 sin 26 , i3n (7 ) for different Arn^j values for 
L = 8000 km as a function of E^. For very long baselines, one can see from fig. || that 
increasing Am^ mainly shifts the resonance energy to higher values thus demanding higher 
beam energies to reach optimal statistical significance. The local maximum in fig. |8] shifts 
to 10 GeV for minimal Am| x and to 40 GeV for maximal Am^. Improved knowledge of 
Amjj would thus in principle allow to discuss optimization issues, but such a study should 
also include systematics and backgrounds. Moreover, the energy distribution of wrong sign 
muon events would add important information to the simple counting of events. A detailed 
discussion of muon energy optimization depends therefore on the way that information will 
be exploited. 



5 Subleading Ara^-effects 



The results shown so far were obtained in the limit Am^ = which is, as already explained, 
a perfect approximation for Am^ <C 10~ 4 eV 2 and/or sin 2 29\ 2 <C 1, i.e. for the VO and 
SMA MSW solutions of the solar neutrino problem. The LMA MSW solution, that we 
will consider in this section, allows however Am^ values up to 2 • 10 _4 eV 2 and prefers 
sin 2 29 12 — 0.8 |J7|, so that effects associated to Am^ can become important, especially for 
the largest Am^ values in the LMA range. Two more mixing parameters, namely 9i 2 and 
5, become relevant if Am 21 is non-negligible [32, 38] (see also, e.g., [39|), as the expression 



3ui 



eq. fljj), including the leading order Am^- corrections shows. While 9 
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for P\ 

is rather constrained, so that we will use sin 2 2# 12 = 0.8 in the following numerical results, 
any value of S in its range < 5 < 2tt is allowed at present. As it is clear from eq. (|16"D, in 
order to calculate the effects associated with a non- negligible Am^, a value of 5 must be 



specified. In fig. [10], the total rates in the appearance channels v e — > (left) and v e 
(right) for Am\ x = (solid line) are compared with the total rates for Am^ = 10 _4 eV 2 
and four possible values of 5 in its range, 5 = 0, tt/2, it, 3tt/2 (dashed lines). As it follows 
from eq. (|l~6| ) and these figures, the size of the effects depends crucially on the value of #13. 
Fig. [10] assumes a value of 6*13 at its upper limit, i.e. sin 2 2^13 = 0.1, whereas fig. [Ll] shows 
the effects for a sin 2 2#i3 one order of magnitude smaller, sin 2 2#i3 
obtained for N.N^e = 10 21 . 



0.01. Both figures were 



Figs. and |TT] illustrate several interesting features of the Am^i effects. First of all, a 
comparison of figs. [It] and [II] confirms that the relative size of the effects grows when 
#13 gets smaller [|J. This is because the zeroth order approximations (in Am^) for the 
appearance probabilities have a sin 2 2#i3 suppression, whereas the linear Am^ corrections 
(CP-conserving and violating) are suppressed by only one power of sin 26*13 (see eqs. (|T0l) and 
(|i~6|) ) and the corrections quadratic in Am^ (CP-conserving) are not suppressed by #13 at 
all. Indeed, at sufficiently small L we have sin(A^ m L/2) = AE m L/2, sm(R + AE m L/2) = 
{R + AE m L/2), etc., and as can be shown p4| , the leading CP-conserving and CP-violating 
Am^- corrections containing the factor cos 9' 12 in eq. (16) reduce essentially to their vacuum 
oscillation form: 
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Figure 10: Appearance event rates n At +(/i~), n fl -(fi + ) in matter with subleading Am^ = 
10 _4 eV 2 and four possible values of the CP-phase 5 = 0, ir/2, n, 3n/2 against baseline L 
(dashed lines) compared with the corresponding event rates with negligible Am^ (solid lines) 
for the channels v e — > z/ M (left) and v e — > (right). Both figures correspond to iV M = 2-10 20 , 
e = 50%, = 20 GeV, sin 2 2# 23 = I, sin 2 26 12 = 0.8 and sin 2 26 13 = 0.1. 
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Figure 11: Same as fig. [JJj but for sin 2 26*i 3 = 0.01. 

where 

Jq P = - sin 5 C13 sin 2#i 2 sin 28 13 sin 26* 23 (37) 
8 

is the CP- violation rephasing invariant of the lepton mixing matrix . 

Unlike the appearance channels, the disappearance channels are dominated by the transi- 
tions to v T {v T ) and are therefore not suppressed by #i 3 or Am^, so that the 6*i 3 -suppressed 
Am^ corrections are much less significant in this case. Figs. [Tl] and |TT] show also that the 
range of the Am 21 corrections is essentially determined by #13, but the precise size and sign 
of these corrections within this range is unknown if S is unconstraint. This can be seen also 
explicitely in figs, pi] and [TT] where the Am 21 corrections show up as oscillations around 
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the leading Am| x contribution to the rates (dashed line), whose initial phase depends on 5. 
The Am| x effects represent consequently in the present LMA scenario in a high statistics 
long-but-not-too-long baseline measurement of #13 an important source of systematic error, 
unless 5 is measured f|, [7|]. On the other hand, sin 5 could be measured or constrained in this 
scenario by comparing the rates in the two CP-conjugated channels u e — » and V e —>■ if 
very high intensity sources and large detectors will become available 0. By comparing e.g. 
the left and right plots of in fig. [TTJ one can see that, as it follows from eq. flTE|), the Am| x 
correction has the same sign in the two channels for the two CP-conserving values 5 = 0, it, 
when sin 5 = 0, but has opposite sign for the two CP-violating values 5 = ir/2, 3n/2, i.e. 
when sin 5 = ± 1 . 

Note, however, that a measurement of 5 based on a comparison of CP-conjugated rates at 
a single value of L could not be enough in order to keep the Am^ effects under control. 
Suppose, for example, that 5 = or 5 = it. Then the Am^ effects would be the same in 
the two channels and an ideal experiment looking for CP- violation would measure sin 5 = 0. 
The value of 5 would, however, still be undetermined for the simple reason that both 5 = 
and 5 = re give sin 5 = 0. As a consequence, it would be still unknown whether the 
Am^ corrections in both channels add (5 = 0) or subtract (5 = ir) to the leading Am^ 
contribution to the rates and this ambiguity would translate, e.g., in an uncertainty on a 
measurement of #13. Such an ambiguity could be resolved by comparing rates measured 



at different distances L. Fig. [TT| for sin 29i 3 = 0.01 shows in the v e — > z/„ channel that 



the change in rates between L = 3000 km and L = 700 km allows to discriminate between 



5 = and 5 = ir. Fig. [10] shows that the different L dependence is also significant for 
sin 2 2013 = 0.1, allowing also to distinguish between the two possibilities 5 = and 5 = n. 
If CP-violation were maximal, | sin<5| = 1, the comparison of the rates at different baselines 
would be less significant but still helpful. Figs. |TU] and |TI] show finally also that the Am^ 
effects become smaller when L is increased. 



6 Matter Effects in the Energy Spectrum. 

Motivated by the small differential event rates, we discussed so far only the influence of 
matter effects on the total rates of wrong sign muon events. We demonstrated in the 
previous chapters that statistical significant deviations from the total event rates in vacuum 
represent already a good test of the MSW theory. A significant test would however be 
also given by a detailed measurement of MSW effects in the neutrino energy spectrum, 
which is modified in a very characteristic way. The —>■ z/ M and — > disappearance 
channels are again dominated by transitions to v T and v T while v e and v e transitions are 
only small corrections. These channels are therefore mostly insensitive to matter effects in 
the differential event rate spectrum and will not be discussed further. 

To understand the effects in the appearance spectrum of v e — > and v e — > we us again 
the approximation Am^ = 0. Matter effects have no influence on the angle #23 in this case, 
whereas they modify the mixing due to #13 significantly. One obtains thus for the v e — > 
and v e —> appearance channels the usual two flavour picture where sin 2# 13 — > sin 28^ = 
sm28 13 /C±, and where C + corresponds to neutrinos, C_ to antineutrinos. The enhancement 
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Figure 12: Modifications in the differential event rate spectrum (events per GeVj due to 
matter effects for = 20 GeV. The solid lines correspond to oscillations taking place in 
the earth mantle, while the dashed lines are for oscillations in vacuum, with Am^ = in 
both cases. The "enhancement" , the "broadening" and the "shift" of the first oscillation 
maximum in the v e — > channel as E decreases, caused by the MSW effect, is clearly seen. 
The dotted lines show for comparison an example o/Am^ corrections for Am^ = 10~ 4 eV 2 ; 
5 = and sin 2 29 12 = 0.8. The assumed parameters are in all cases N^N kT e = 10 21 , 
L = 6596 km (CERN-MINOS), Am 2 31 = 3.5xl0~ 3 eV 2 and sin 2 26 13 = 0.01. 



of 8'^ is maximal in the neutrino channel when the neutrino energy E coincides with the 
MSW resonance energy E ies defined via: 



ff-os^. (38) 



Note however that the maximum of the event spectrum in the u e — ► channel does in 
general not coincide with E res since the MSW oscillation probabilities are folded with the 
fluxes and cross sections. The maximum of the event spectrum is thus determined by the 
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Figure 13: Same as in fig. [7^ but for = 50 GeV. For details see text. 
maximization of 

p^)) 2 W£/£,), (39) 

but the resulting maximum is for the muon energies under consideration still around E res . 
The offset depends in a rough approximation on the difference between E Tes and the maxi- 
mum of the flux which lies roughly at an energy of the order E^. 

This has to be compared with the vacuum case where the oscillation probabilities are also 
folded with fluxes and cross sections and where the maxima of the event rates are also not 
precisely at the maxima of the oscillation probabilities. The event rate spectrum is thus 
due to this folding in both cases with and without matter a rather complicated function of 
E v which depends in a non-trivial way on L, E^ and Am^. Nevertheless for given L, E^ 
and for given Am^, 6*23, #13 measured with a suitable long baseline experiment, one can 
predict the shape of the differential event rate spectrum in all channels and compare it with 
the spectrum of oscillations unaffected by matter. This results in a very good opportunity 
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to detect specific details of the MSW effect which arise when the oscillation parameters 
are chosen such that the first maximum of vacuum oscillation coincides roughly with E res . 
The point is that the MSW effect changes the probabilities compared to vacuum in three 
genuine ways: The first maximum of the oscillation probability as the energy decreases is 
enhanced, its width is broadened and its center is shifted to lower energies. Similarly one 
has an "anti-MSW effect" in the antineutrino appearance channel which implies for the first 
oscillation maximum a reduction in height, again a broadening and a shift to lower energies. 

These "genuine MSW effects" are demonstrated in fig. [L2| (where = 20 GeV) and fig. [13| 
(with = 50 GeV) showing the modifications in the energy spectrum (events per GeV) 
due to matter effects. The solid lines are in matter while the dashed lines are without matter 
and the assumed parameters are NpNvre = 10 21 as before, L = 6596 km and sin 2 26 13 = 
0.01. Fig. |I2| shows already all effects due to the MSW mechanism, namely the broadening 
(the last oscillation in matter covers two oscillations in vacuum), the shift (the maximum 
in matter lies almost in the minimum in vacuum) and the enhancement or suppression 
compared to vacuum. Figs. [12| and [13| also include an example of Am^ corrections to the 
spectrum in the LMA scenario (dotted lines). Size and sign of the corrections depend on 
the value of the phase 5. The dotted line assumes 5 = and moreover Am^ = 10 _4 eV 2 
and sin 2 29^ = 0.8. As already observed in the previous section, the Am^ corrections are 
relatively small for the very long baselines considered here. 

It is interesting to look at the modifications when the muon energy becomes higher, e.g. for 
Efj, = 50 GeV as shown in fig. |H| The point is that the beam energy is already rather far 
away from the MSW resonance energy and the importance of the weight function f„ eV (and 
of the scaling of unoscillated events) in the determination of the shape of the spectrum 
increases. Thus the genuine broadening, shift and enhancement / suppression effects become 
harder to distinguish. For the u e — > channel the effect could be hard to distinguish from 
uncertainties (with low statistics) in the spectrum. This brings up the general issue that one 
has to have enough statistics for such an analysis. In order to have a chance to see such effects 
one has to be lucky and 6*13 should be at the upper experimental limit (see scaling laws). 
Otherwise N^N^e must be increased correspondingly which implies a more intense muon 
source, a larger detector or both. Although a more quantitative analysis of the significance 
of effects in the differential neutrino event rate spectrum is beyond the scope of this paper, 
an analysis of the differential event rate spectrum would clearly provide extremely valuable 
additional information which would allow to test some of the characteristic features of the 
MSW mechanism. 



7 Conclusions 

Assuming three-neutrino mixing we studied in this paper the possibility to test the MSW 
effect in terrestrial very long baseline neutrino oscillation experiments which become possible 
with neutrino factories. Such direct tests are important since the MSW mechanism is widely 
used in different scenarios of neutrino physics and astrophysics. The correct analysis and 
interpretation of the data from terrestrial very long baseline neutrino oscillation experiments, 
L ^ 1000 km, is in fact impossible without a proper treatment of matter effects. The latter is 
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also crucial for the searches of CP-violation in neutrino oscillations generated by the lepton 
mixing matrix, since matter effects create an asymmetry between the two CP-conjugated 
appearance channels. 

Studies of the v e —>■ and v e — > oscillations are by far most promising for the detection 
of the matter effects since the corresponding total event rates are affected in a drastic way 
by these effects. We considered for the present study neutrino trajectories through the 
earth which cross the mantle, but do not pass through the earth core, which corresponds 
to neutrino path lengths L ^ 10600 km. Using analytic expressions for the three neutrino 
oscillation probabilities in matter in the constant average density and small Am| approx- 
imations and including fluxes, cross sections and detection efficiencies allows to describe 
the relevant event rates analytically. This permitted a full analytic understanding of our 
numerical results. 

By considering the asymmetry between the v e — > and v e — > induced wrong sign muon 
event rates, we studied the statistical significance of the observation of matter effects as a 
function of the neutrino oscillation parameters #13 and Am| x as well as its dependence on 
the experimental conditions via the parent muon beam energy E^, the path length L and 
the product of useful muons, detector size and efficiency The scaling of rates, statistical 
significances and sensitivities with the relevant mixing angles, in particular, with #13, the 
intensity of the muon source and with the detector size and efficiency have been given, 
so that the results for any value of those parameters can easily be obtained. The sign of 
the asymmetry depends on whether the two closest neutrino mass eigenstates are lighter 
(Atojj > 0) or heavier (Am^ < 0) than the third one, thus providing a way of determining 
which of these two possibilities is realized. Figs. |6| and |7| show the conservative ranges of 
sin 2 2#i3 where that determination would be significant at a given confidence level from the 
statistical point of view as function of the baseline. 

We analyzed, in particular, the statistical significance of matter effects as a function of 
and L. The most important requirement regarding the muon energy is that for given Am| x 
it has to be greater than the MSW resonance energy. For, e.g., Am^ < 6.0 (8.0) x 10~ 3 eV 2 , 
this implies E^ ^ 20 (30) GeV. Using a simple "total rates based" method and assuming 
that #13 is known, any L from the interval (4.0 — 10.0) x 10 3 km gives relatively good 
sensitivity to matter effects if Am 2 31 = 3.5 x 10" 3 eV 2 , sin 2 29 r3 = 0.01 and E^ = 20 GeV, 
with the sensitivity increasing with L. For = 50 GeV, the sensitivity varies very little 
for L = (4.0 - 10.0) x 10 3 km. Let us note that if Am 2 31 ~ 10~ 3 eV 2 and sin 2 26 n £ 0.01, 
then establishing the matter effects (or obtaining a stringent upper limit on sin 2 2^13) would 
require ~ (40 — 50) GeV. As our results show, the optimal value of depends on the 
precise value of Am 31 . 

Our analysis shows that a higher sensitivity to the MSW effect in the case of relatively 
small values of sin 2 2# 13 is achieved at very large L. We showed that this conclusion holds 
also when subleading Am|i effects are included. These effects can be significant in the 
case of the LMA MSW solution with Am\ x ~ (0.5 - 2.0) x 10" 4 eV 2 . At very large L, 
the indicated Am^-induced effects are essentially negligible and the results obtained in 
the limit of Am^ = are therefore sufficiently accurate. This is valid even in the case 
of the LMA MSW solution with Am 2 x ~ (0.5 - 2.0) x 10~ 4 eV 2 . For shorter baselines, 
Am 2 x effects are non-negligible and a determination of the CP-violating phase 5 would be 
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necessary in order to know their precise magnitude. We have found that the L dependence 
of the Amji effects offers the possibility to determine the CP-phase S, especially when sin 5 
is small. For sin 2 2# 13 = 0.01, for instance, the event rates due to the v e — > and v e — > 
transitions change considerably when L changes from ~ 700 km to ~ 3000 km. Thus, a 
measurement of these rates, e.g., at the indicated two distances could allow to determine the 
value of 5 with a certain precision. Finally we discussed the matter effects in the differential 
event rate spectrum as a function of the neutrino energy, and showed that they lead to very 
characteristic distortions. The observation of these distortions would allow very detailed 
tests of the MSW theory. 

To conclude, our study shows that the predictions of the MSW theory can be tested in a 
statistically reliable way in a large region of the corresponding parameter space by a simple 
analysis of the total event rates in a very long baseline neutrino oscillation experiment. This 
can allow to determine the sign of Am^ as well, as was recently noticed also in ||. Not 
seeing the matter effects would lead to impressive upper limits on the mixing angle #13 down 
to sin 2 26*13 — 10~ 4 or even better. 
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